A theoretical analysis of the chemical bonding and electronic structure of graphene 
interacting with Group lA and Group VIIA elements 
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We propose a new class of materials, which can be viewed as graphene derivatives involving Group 
lA or Group VIIA elements, forming what we refer to as graphXene. We show that in several cases 
large band gaps can be found to open up, whereas in other cases a semimetallic behavior is found. 
Formation energies indicate that under ambient conditions, sp^ and mixed sp^/sp^ systems will 
form. The results presented allow us to propose that by careful tuning of the relative concentration 
of the adsorbed atoms, it should be possible to tune the band gap of graphXene to take any value 
between and 6.4 eV. 

PACS numbers: 71.15.-m, 71.15.Mb, 71.20.Nr 



A new derivative of graphene was recently reported[l|, 
where hydrogenation turned graphene, a recently discov- 
ered two-dimensional C based material, 0-0 into what 
is referred to as graphane. This compound was actu- 
ally studied theoretically before the material was synthe- 
sized, and the materials phase stability as well as elec- 
tronic properties (in particular the presence of a gap) 
were predicted from first principles theory 0, Q. In 
the GGA-based calculations of Ref0, H it was found 
that graphane has every C atom attached to a H atom, 
and that the H atoms arc bonded alternating on both 
sides of the C layer. Furthermore, it was found that 
graphane is a semiconductor with a rather wide energy 
gap (3.5 eV). Subsequent calculations based on the GW 
approximation, 0, which are supposed to result in 
more accurate band gaps, reported on a gap- value of 5.2 
eV.dll 

In the experimental study [l[ graphane was indeed 
found to have a band gap of the electron states. Hence, 
it was demonstrated that the adsorption of hydro- 
gen turned highly conductive graphene into insulating 
graphane. The reason a band gap opens up when hy- 
drogen is adsorbed on graphene is that sp^ bonded C 
atoms can achieve a higher degree of sp'^ bonding. In this 
case three of the four covalent sp'^ bonds are saturated 
by C atoms and the fourth covalent bond is saturated 
by an H atom. The reason for the change in conduc- 
tivity is in accordance to expectations for sp'^ bonded 
carbon (e.g. diamond), which is an insulator. On an 
electronic structure level the transition to an insulator 
is connected to the fact that the p^ orbitals, which in 
graphene form conducting states (called tt and tt* ) at the 
Fermi level, participate in graphane in strong covalent a 
bonds, with a gap between bonding and anti-bonding 
states. This finding is illustrated in the theoretical work 



of RefsIH, and from the adsorption of single hydrogen 
atoms on graphene. [l^ Moreover, in Ref. [13, the scat- 
tering of electrons in graphene by clusters of impurities 
is studied and for reviews on experimental work on alkali 
metals on graphite, see Refs. Mnn 

The possibility to open up a band gap in graphene 
derivatives, like graphane, is very important when consid- 
ering the potential of this material for electronics appli- 
cations. A natural question which then arises is if other 
atomic species of the Group lA elements can produce a 
similar behavior, potentially with a slightly smaller band 
gap, which is more suitable for electronics applications. 
It is also possible that Group VIIA atoms can adsorb 
on graphene, to form sp"^ bonds with the C atoms and 
with a gap in the electronic structure. For the case of 
single impurities, the difference in chemical bonding has 



been analyzed in Ref. [18|. It was found that alkali metals 



(except Li) form purely ionic bonds, with very small dif- 
ference in total energies between different positions (on 
the top of a carbon atom, in the middle of the bridge 
between two C atoms, in the center of a C-hexagon), and 
the same for halogens (except F). This seems to be in 
agreement with the present results on complete cover- 
age. We have here investigated the possibility of a gaped 
electronic structure by performing ab-initio calculations 
of the geometry and electr onic structure for the Group 
I A and Group VIIA elements adsorbed on graphene, and 
we will refer to this new material as graphXene (where 
X represents a Group lA or Group VIIA element). The 
theoretical calculations have been done using the gener- 
alized gradient approximation (GGA) and the GW ap- 
proximation. For simplicity we have only considered the 
so called chair geometry [Tj in our calculations, i.e. with 
one adatom attached on each C atom, alternating being 
above or below the C atoms (for a figure of the geome- 
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try see Fig.l of RefQ)- We have considered several cases, 
for instance with Group lA elements on both sides of 
the C network, and with Group VIIA elements on both 
sides of the C network. We have also considered mixed 
cases where Group lA atoms are above the C network 
and Group VIIA atoms are below the C network. 




FIG. 1: (Color online) Geometry of the C - X bonding (X is 
a Group lA or Group VIIA element, yellow atom, C atoms 
are blue). The bond-angles for the C-C and C-X bond are 
defined as (f>c-c and <j)c-x, respectively. 



We have used the VASP[l9| (Vienna Ab-initio simu- 
lation package), implementing the projector augmented 
waves (PAW) method [20|. to compute the ground state 
geometry and the electronic structure. We obtained a 
reliable geometry of the structure by optimizing atomic 
positions, using the generalized gradocnt approximation 
(GGA). For all the calculations, we have used a 500 eV 
cut-off for the wavefunction. During the optimization of 
the structures, a k-point grid[2l| of 12 x 12 x 1 was used. 

Density functional theory (DFT) is a ground-state the- 
ory and can therefore not treat excited states. [22j 
Equally important is that DFT does not differentiate 
between the two potentials Vjv-i and Vat that should 
be seen by occupied- and unoccupied states, respec- 
tively. Moreover, self-interaction is not included and 
Yxc does not include non-local effects nor energy- and 
electron density dependencies. To calculate accurate 
band gaps a more reliable method is needed and here 
we use the all-electron GW implementation based on the 
FP-LMTO method by ReflH. The GW approximation 
of Hedinji, [l^ can be derived in a systematic way from 
many-body perturbation theory, and it is common to take 
an LDA or GGA calculation as the starting point. The 
result from the theoretically complex derivation is 
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which is the quasi-particle energy resulting from first 
order perturbation expansion of an operator, that is 



made from the difference between the self-energy oper- 
ator (E = iGW) and the LDA exchange-correlation po- 
tential (V^c^)- G and W is the LDA Green's function 
and screened Coulomb interaction, respectively. W is re- 
lated to RPA polarization via the dielectric function. In 
this work a one-shot GW is used and we have 
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The quasi-particle renormalization factor Zi^k is given by 



l-(z,fc| 



and is around 0.9 in our calculations. 

There are several reviews of the GWA method in the 



literature |2J-|27[. Concensus is that GWA improve the 
bandgap of semiconductors and insulators and the calcu- 
lated GWA bandgap is often within 10% of experiment. 

In general one may expect that depending on the el- 
ement which is attached to the C atom of the graphene 
network, a varying degree of sp'^ and sp^ bonding will oc- 
cur. The best way to characterize the degree of sp^ or sp^ 
bonding may be to use the bond angles between C-C of 
the graphene network as well as the bond angle between 
the adatom and the C atoms, i.e. the C-X angle. These 
angles are shown in Fig.l, where the C-C angle is labeled 
(f>c-c and the C-X angled is labeled (pc-x- For a purely 
sp^ bonded system 4>C-c is 120 degrees and (pc-x is 90 
degrees. For a purely sp"^ bonded system (pc-c is the 
same as (pc-x, with a value of ^ 109.5 degrees. For in- 
termediate cases one will observe bond angles somewhere 
in between these two extreme cases. 

In Table I we list the calculated bond angles for all the 
systems investigated in this study. It may be seen that 
elemental adsorption with the same chemical species on 
both sides of the graphene sheet produce a situation with 
ideal sp^ bonding for CI, Br and I adsorption. In the 
cases of adsorption of Na, K, Kb and Cs one observes a 
chemical bonding which is almost of pure sp^, whereas 
for Li one observes an even mixture between sp^ and sp^ 
bonding. When H or F are adsorbed, bonds which are 
mostly of sp^ character are observed. Among the mixed 
systems, the combination of H on one side and F on the 
other side, or CI on the other side, produce a situation 
with dominant sp'^ bonding. We will below refer to these 
cases as HF and HCl adsorption, respectively. A combi- 
nation of Li and F, as well as Li and H, results in similar 
scenarios, but with a higher degree of sp^ bonding. It is 
interesting to note that Li and H, and possibly Li and 
F appear to be semi-metals. Judging from Table I, we 
speculate that the transition in bond angle from semicon- 
ducting (or insulating) to semi- metallic is at around 114° 
for the C-C-C bond. Finally we note that a combination 
of H and Br produce a situation which reflects a bonding 
closer to pure sp^ character. 
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C2H2 


2.54 


1.54 


1.11 




111 


107 


sp^ 


3.49 


5.74 


-0.11 


C2Li2 


2.58 


1.53 


2.02 




116 


77 


mix 


wm 




-0.64 


C2Na2 


2.67 


1.54 


2.63 




120 


88 


sp2 


m 




0.08 


C2K2 


2.82 


1.65 


2.67 




118 


81 


sp2 


m 




0.97 


C2Rb2 


2.86 


1.67 


3.01 




118 


81 


sp2 


m 




1.29 


C2CS2 


2.94 


1.72 


2.92 




118 


82 


sp2 


m 




1.79 


C2F2 


2.61 


1.58 


1.38 




111 


108 


sp3 


3.10 


7.4 


-0.81 


C2CI2 


2.57 


1.47 


3.81 




121 


90 


sp2 


m 




0.41 


C2Br2 


2.69 


1.55 


4.18 




120 


90 




m 




0.58 


C2I2 


2.87 


1.66 


3.99 




120 


90 


sp2 


m 




1.21 


C2HF 


2.57 


1.56 


1.10 


1.39 


111 


108 


sp^ 


3.11 


6.38 


-0.47 


C2HCI 


2.77 


1.66 


1.10 


1.76 


112 


107 


SP^ 


0.87 


2.91 


0.41 


C2HX" 


2.65 


1.60 


1.39 


1.80 


112 


107 


SP^ 


0.26 




-0.08 


C2HBr 


2.68 


1.57 


1.13 


3.70 


118 


98 


sp^ 


m 




0.72 


C2LiF 


2.59 


1.55 


2.16 


1.44 


114 


105 


mix 


wm 




-0.39 


C2LiH 


2.58 


1.53 


2.19 


1.14 


115 


104 


mix 


sm 


0.30'' 


0.01 



''X=C1, F alloyed 50-50 

''Using GGA gives semi-metallic behavior for C2LiH and with the 
GW correction a small in-direct gap opens up between K (VB) and 
r (CB). 

TABLE I: Geometrical and electronic structure data for all 
systems examined. Bond distances are given in [A], angles 
in [°], lattice constant (a) in [A], gap energies in [eV], and 
formation energy in eV/atom. Interatomic bond lengths are 
denoted C-C between carbon atoms and C-X (or C-Y) for 
other bond distances. The systems that show a band gap 
are all direct gap. C2LiH is a semi-metal (sm), C2LiF is a 
borderline case sm/weak-metal and rest are metallic (m) or 
weak-metals (wm) at the GGA level. 



The energy of formation (E''^ in Table I) is cal- 
culated using E^=Et(graphXene)-Et(graphene)-Ej(X)- 
Et(Y) where Ej is the total energy, graphXene=C2XY 
where X,Y is one of H, Li, K, Na, Rb, Cs, F, CI, Br or 
I. Note that Et(graphene) is the total energy per unit 
cell, i.e. for two C atoms. It is seen from Table I that 
all sp^ systems except C2HCI have negative energies of 
formation, indicating that the compounds will form. The 
material C2HCI will form if alloyed with 50% F (lower F 
concentrations have not been investigated). 

Moreover, if the reference state of the X or Y atoms 
is not the gas phase at ambient conditions, but a gas 
phase at elevated pressure (P) and temperature (T), a 
term proportional to k^T/jiP, enters the expression for 
the free energy, in such as way as to stabilize the phase 
with X and Y adsorbed on the graphene layer [2^. Hence 
it is possible to stabilize some of the materials listed in 
Table I, even if E'' > 0. 

In the cases with dominant sp"^ bonding it is expected 
that a gap opens up, whereas for the intermediate cases 
the situation is less clear. The electronic structures of the 
presently studied systems show that a band gap opens up 
for H adsorption on both sides of the carbon layer as well 
as for F adsorption on both sides. The band structure 
of H adsorption has been shown before Q and for this 
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FIG. 2: (Color online) Band structure of graphene with F 
on both sides of the C layer, in a chair geometry using GGA. 
Circles show the GW corrected energies. The rigid shift of the 
unoccupied states is seen when the lowest GGA conduction 
band is plotted but shifted to fit the GW values (dotted line 
with crosses). 
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FIG. 3: (Color online) Band structure of graphene with H 
on one side and F on the other of the C layer, in a chair 
geometry, using GGA theory. Circles show the GW corrected 
energies. The rigid shift of the unoccupied states is seen when 
the lowest GGA conduction band is plotted but shifted to fit 
the GW values (dotted line with crosses). 



reason we do not repeat this electronic structure here. 
However, we show the calculated energy bands of F ad- 
sorption in Fig. 2. The general shape of the energy band 
structure is the same for the GGA and GW calculation, 
except for a nearly rigid shift of the unoccupied states 
in the GW calculation, so as to open up a larger band 
gap. This nearly rigid shift is illustrated in Fig. 2 (for F 
on both sides of the graphene layer) and Fig. 3 (for F on 
one side and H on the other side of the graphene layer), 
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by comparing quasi-particlc GW-energies (circles) with 
the energies obtained by shifting rigidly the bottom of 
the conduction band from the GGA calculation (crosses) . 
Both materials show a shifted GGA band structure which 
coinsides with the GW bands. This shift is in accordance 
with Fig. 1 (for diamond) of the ground-breaking paper 
by Hybertsen and Louie [2^. An additional evidence 
for the nearly rigid shift in the GW calculation is illus- 
trated in Fig. 4, which shows the quasi-particle energies 
plotted against the Kohn-Sham eigenvalues for C2HF. 
From least-squares fitting (full lines) the valence bands 
and conduction bands are found to have the slopes 1.00 
and 1.14, respectively. A slope of 1.0, with no scatter of 
the data around the fitted line, would indicate a perfect 
rigid shift for electron states of any k-point, and the data 
in Fig. 4 gives good evidence for that the GW states are 
shifted more or less rigidly. 
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Kohn-Sham eigenvalue [eV] 



FIG. 4: (Color online) Quasi particle energies plotted against 
the Kohn-Sham eigenvalues for C2HF. The full lines show the 
least squares fit for the valence bands (VB) and conduction 
bands (CB), respectively. 



non trivial, but on the other hand, this leads to more 
flexibility to obtain the desired value for the band gap, 
a key quantity in view of using graphene in electronic 
devices. 

In conclusion, we have demonstrated, using accurate 
GGA and GW calculations, that a reasonable band gap 
can be opened up for graphene derivatives with Group 
I A and Group VII A elements, forming a new class of 
materials we refer to as graphXene. Depending on ad- 
atoms, graphXene forms sp^ bonds, sp'^ bonds or a mixed 
sp^/sp'^ binding. Several of the proposed materials have 
formation energies which suggest a stable configuration, 
where in general the sp'^ and mixed sp^/sp'^ systems 
should be experimentally observed. We point out that 
although not all systems investigated here have a forma- 
tion energy at ambient conditions which suggests that 
they should form, it is likely that elevated temperature 
and pressure will stabilize them. 

The most important conclusion of our study is that a 
careful mixing of different dopants of the graphene layer 
enables a tuning of the band gap from to ^ 6.4 eV, a 
result which has obvious technological ramifications. For 
the sp'^ bonded, gaped systems, the C-C-C angle is close 
to 109.5°. For the sp^ bonded systems the G-G-G angle is 
120°; or just below this value. The intermediate cases are 
semi-metallic and we speculate that the semiconducting 
(insulating) semi-metallic transition appears for a C- 
C-C angle which is around 114°. 
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For the adsorption with different chemical species on 
the two sides of the C layer, both H-C-F and H-C-Cl are 
found to display band gaps. Since the two materials have 
electronic structures which are rather similar, we show 
only the GGA electronic structure of H-C-F in Fig. 3. The 
H-C-F system has a direct, 3.11 eV band gap, as is seen 
from Fig. 3. In GW theory this band gap is 6.38 eV. The 
H-C-Cl case results in a smaller direct, band GGA gap 
of 0.87 eV, which opens up to 2.91 eV from GW theory. 
Therefore, a whole range of bandgaps can be induced in 
graphene by carefully choosing the chemical species of 
the adatoms. One can even go further by using not only 
two kind of dopants, but three: the alloyed structure H- 
G-(Gli/2Fi/2), gives a GGA gap of 0.26 eV, which is very 
different from 0.87 eV of H-C-Cl and 3.11 eV of H-C-F. 
This result demonstrates that the relationship between 
the different elements concentration and the bandgap is 
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